Microscopic Viscoelasticity: Shear Moduli of Soft Materials 
Determined from Thermal Fluctuations 
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We describe a high-resolution, high-bandwidth technique for determining the local viscoelasticity 
of soft materials such as polymer gels. Loss and storage shear moduli are determined from the 
power spectra of thermal fluctuations of embedded micron-sized probe particles, observed with an 
interferometric microscope. This provides a passive, small-amplitude measurement of rheological 
properties over a much broader frequency range than previously accessible to microrheology. We 
study both F-actin biopolymer solutions and polyacrylamide (PAAm) gels, as model semiflexible 
and flexible systems, respectively. We observe high-frequency oj'^^* scaling of the shear modulus in 
F-actin solutions, in contrast to lo^^^ scaling for PAAm. 

PACS numbers: 61.30.Cz, 61.30.-v, 64.70.Md, 83.70.Jr 



The local analysis of viscoelasticity can explore the 
small-scale structure of complex fluids. It can also po- 
tentially measure bulk viscoelastic quantities in small 
samples. Experiments have been done to characterize 
the local viscoelasticity of materials since at least the 
1920s, when magnetic particles in gelatin were manipu- 
lated by field gradients . Crick used a similar technique 
to study living cells |^ . With the recent advent of meth- 
ods for force generation, detection, and manipulation of 
particles on sub-micrometer scales, experimental possi- 
bilities have expanded greatly, and interest in microrhe- 
ology has grown substantially |^^]. Here we report an 
optical technique for high resolution and high bandwidth 
observations of the thermal fluctuations of particles em- 
bedded in soft materials. Using dispersion relations from 
linear response theory, the frequency-dependent loss and 
elastic storage shear moduli are both calculated from the 
fluctuation power spectra, which we measure with a res- 
olution of 2 Afrom 0.1 Hz to 20 kHz. This frequency 
range exceeds that of video-based microrheological ex- 
periments (although large-volume rheometers can reach 
frequencies up to 500 kHz [^) and allowed us to observe 
high-frequency ur"/^ scaling of the shear modulus for en- 
tangled, semiflexible solutions. These dynamics differ 
fundamentally from those of flexible polymer systems. 

We have studied F-actin solutions as a model semi- 
flexible polymer, and polyacrylamide (PAAm) gels as a 
flexible polymer control. Actin is one of the primary 
components of the cytoskeleton of plant and animal cells 
and is largely responsible for the viscoelastic response of 
cells ||l^ . Actin is a particularly accessible model system 
because individual fllaments can be hundreds of microns 
in length. Viscoelastic properties of entangled F-actin 
solutions in vitro have been measured using conventional 
macroscopic rheology piHl5[. Actin has also been the 



subject of recent microrheological studies 

The relationship between thermal motion and hydro- 
dynamic response is well-known in the context of Brow- 
nian motion in simple viscous fluids. Less obviously, the 
bulk viscoelastic properties of complex fluids can also be 
determined from thermal motion. Such a passive tech- 
nique is particularly useful for biopolymer systems which, 
on the one hand, have a limited range of linear response, 
but which, on the other hand, have a small enough shear 
modulus that the thermal motion of micron-sized probes 
is detectable. This idea was exploited recently by Mason 
and Weitz iQ, who used diffusing wave spectroscopy 
(DWS) to observe average thermal fluctuations of con- 
centrated suspensions of probe particles and to measure 
rheological quantities in complex fluids to high frequen- 
cies. That technique measures light multiply-scattered 
by many embedded particles, assumed to be independent. 
While DWS requires large samples, the present technique 
monitors isolated particles in very small samples. 

For spherical probes, the response is given by a gener- 
alized Stokes formula, which allows for quantitative de- 
termination of macroscopic viscoelastic moduli. In our 
experiments, rigid spherical beads of radius R move in a 
viscoelastic medium consisting of polymer plus solvent. 
When a force / acts on such a sphere, the deformation of 
the surrounding medium can be calculated ||l7|| to obtain 
the linear response for the sphere: = a{u})f^, where 
X is the position of the sphere. In general, the response 
function is complex: a(w) — a'{uj) -\- ia"{uj). Within 
certain limits described below, this response function re- 
flects pure shear motion of the medium and is given by 
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Here, the frequency-dependent complex shear modulus is 
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given by G{uj) = G'{lu) + iG"{u}), where G' and G" are 
the storage and loss modulus. (For the case of a purely 
viscous fluid G = —iurj, leading to the well-known Stokes 
formula / = GnriRx). The fluctuation-dissipation theo- 
rem ijl^ relates a, and thus G, to Brownian motions of 
the bead induced by thermal fluctuations of the medium. 

As in the case of a viscous fluid Eq. (|l|) is valid 
provided that the viscous penetration depth S — y^2T]/poj 
is large compared with R, where p and 77 are solvent 
density and viscosity; this is true for frequencies up to 
about 1 MHz for a 1 pm bead in water. The validity of 
continuum viscoelasticity also requires that the bead is 
at least larger than the mesh size of the network. 

A somewhat more subtle limitation on all microrhe- 
ological methods to date, including ours, concerns the 
two-component nature of polymer solutions. In a gen- 
eral viscoelastic medium, there can be another relevant 
material parameter in addition to the shear modulus: 
the compression modulus, or equivalcntly, the Poisson 
ratio. In a simple viscous fluid, only shear stress is rel- 
evant because of the incompressibility of the fluid. In a 
polymer solution consisting of both network and solvent, 
this is also true for the high-frequency response because 
the network and incompressible solvent become strongly 
coupled. While a rigorous calculation of the response 
requires a full treatment including both network and sol- 
vent aspects of the material, we can estimate the validity 
of the approach given above within a so-called two fluid 
model [Q, which takes into account the viscous cou- 
pling of the two components. We estimate the frequency 
crossover from the balance between viscous and elastic 
forces. The viscous force of the solvent on the network 
is rjv/S,'^ per volume, where ^ is the mesh size and v is 
the solvent velocity relative to the mesh. The local elas- 
tic force in the network is GV^u ^ Gu/R^ per volume 
at the bead surface, where u is the network displacement 
field. Thus viscous coupling is dominant, and solvent and 
network move together, above a frequency 



the Kramers-Kronig relations by evaluating the dis- 
persion integral 
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Above ujc, Eq. (|l|) will hold. In our experiments Uc is 
of order 10 Hz (G ~ 1 Pa, ^/R =i 0.1, rj = 1 cP), which 
is at the lower limit of our instrumental resolution. At 
frequencies uj <^ ojc, the response includes compressional 
modes not present in Eq. (|^); these will have wavelengths 
q^^ ^ ^{G/rjLuY^^ and will make a(w) independent of 
bead size. 

Applying the fluctuation-dissipation theorem to the 
power spectral density (PSD) of bead displacements (x^) 
yields 
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When a"{uj) is obtained over a large enough range of 
frequencies, the complex response a{uj) is obtained from 



The principal-value integral P in Eq. (^ is conveniently 
computed on a discrete data set as a sine transform of 
a" {ijj) followed by a cosine transform. We finally take 
the complex reciprocal of a{Lo) to obtain the complex 
G{uj) by Eq. (|l|). We have evaluated Eq. (^) using model 
spectra with various upper and lower cutoffs |171| and 
found it to be reliable over the scaling range discussed 
below. 

The shear modulus G{uj) typically exhibits three dis- 
tinct regimes in entangled flexible polymer solutions . 
For a non-crosslinked solution the behavior is essen- 
tially that of a viscous liquid at frequencies lower than 
where tr is the reptation time. Above I/tu a 
rubber-like plateau appears, with a dominant, frequency- 
independent elastic response as in a crosslinked gel. At 
even higher frequencies, above the relaxation rate of the 
mesh 11^], the high-frequency moduli G' and G" are ex- 
pected to obey a pure power-law, since no new character- 
istic time appears until the molecular cutoff is reached. 

Our samples were prepared as follows. Actin was ex- 
tracted from chicken skeletal muscle according to stan- 
dard procedures p^ , combined with silica beads, and 
polymerized (buffer: 2 mM HEPES (pH 7.0), 2 mM 
MgCl2, 50 mM KCl, 1 mM EGTA, 1 mM ATP) at a 
concentration of 2 mg/ml under slow rotation to avoid 
bead sedimentation. This yields a highly entangled so- 
lution of stiff fllaments with persistence lengths of order 
10 ^m and a mesh size ^ = 0.2 /j,m 23 1. PA Am gels were 
prepared as in standard gel electrophoresis |2^ with con- 
centrations of 2, 2.5, and 3% weight/ volume (20, 25, and 
30 mg/ml); the PAAm contained 3% bis-acrylamide as 
a crosslinker. These flexible polymers have a persistence 
length on the order of nanometers; the mesh size of a 
2% gel is about 50 A |2^]. The smallest sample dimen- 
sion was their thickness of 70 /im (F-actin) and 140 /xm 
(PAAm). 

For bead position detection we couple a near-infrared 
laser into an optical bench-mounted differential interfer- 
ence contrast (DIG) microscope. Polarized at 45° to the 
axis of the DIG WoUaston prism, the beam is split and 
the microscope objective produces two overlapping foci 
(~200 nm apart) of orthogonal polarizations |^^. A 
bead of higher refractive index than the medium in the 
double focus causes a differential phase shift A(^(a;) be- 
tween the beams which is a measure of bead position 
X. After recombination by a second WoUaston prism, 
the light, now elliptically polarized, is passed through 
a quarter-wave plate and a polarizing beam splitter onto 
two photodiodes. The normalized difference of the photo- 
diode currents is sin(A(/!)). Small excursions of a centered 
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bead result in a proportional voltage response which we 
calibrated by controlled displacements of a piezoelectric 
stage. This set-up is also used for optical trapping; here, 
the laser focus was broadened and the power reduced 
(~0.6 mW in the specimen) so that optical forces were 
negligible. Our resolution at high frequencies was limited 
by shot noise to about 2 x 10~^ nm^ Hz~^. The bead- 
position signal was digitized at 60 kHz and filtered above 
20 kHz to eliminate aliasing. 

In actin solutions, the power spectral density (PSD) of 
bead fluctuations (x^) was highly reproducible in shape 
and amplitude, not only for different beads within a 
sample but between samples and between different actin 
preparations |17|. For samples with beads of diameter 
0.9, 2.1, and 5.0 fim, the variation in PSD among six 
beads per sample (at 100 Hz) was 10 — 20%. 

The fluctuation power spectra of beads in actin are 
shown in Fig. |l|. Above about 20 Hz, the PSDs scale 
(roughly) with bead size, as predicted by Eq. (|l|). All 
spectra show a slight downturn above 3 kHz, which is 
currently an experimental limitation. At low frequencies, 
the slopes decrease in a bead-size dependent way. We be- 
lieve this to be the transition from a regime of pure shear 
fluctuations to long-wavelength compressional modes, as 
predicted by Eq. (||). For 5.0 /im beads, the PSD dis- 
plays power-law behavior uj~" with n = 1.77 ± 0.02) 
over almost three decades of frequency. In contrast, 
0.8 /im polystyrene beads in glycerol, a simple viscous 
fluid, should display a diffusional exponent oin = —2; we 
find n — 1.94±0.03. Furthermore, from this spectrum we 
obtain a value of 11.3 P, compared to a tabulated value 
of 9.3 P at 25°C. 

The G'{uj) and G"{uo) calculated from the spectra for 
5.0 /im beads are shown in Fig. |I|B. These beads are 
about 20 times larger than the mesh size, so that the con- 
tinuum model of Eq. (||) applies: as expected, a power- 
law has emerged in G' and G" . No low-frequency plateau 
is observed. (In |jlj] this plateau was found below about 
0.01 Hz at similar concentrations.) The scaling corre- 
sponds, via Eqs. (|]) and (^, to a complex power-law 
G(w) oc {iujY, so that G" /G' = tan7rz/2. This ratio 
allows an estimate of z = n — 1 = 0.76 ± 0.02 between 
10 Hz and 100 Hz, consistent with the power observed in 
the PSD. 

Other studies of actin rheology UJl^, at frequencies 
below about 3 Hz, have reported G ~ lo^/"^. However, 
at these lower frequencies, their results may reflect an 
incomplete transition to the uj'^/^ scaling regime. The 
magnitudes we measure for G' and G" are consistent 
with other experiments &E|,^, although discrepancies 
persist in the literature [|12|. Amblard et al. found 
bead-diffusion dynamics in more dilute actin solutions 
consistent with our scaling results, but suggest an expla- 
nation based on isolated filament dynamics. We believe 
that at our higher concentrations, for which the mesh 
size is substantially smaller than the bead diameter, a 



continuum elastic approach as described here is correct. 
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FIG. 1. (A): Power spectral density (PSD) of thermal bead 
motion in an actin solution (2 mg/ml), averaged over several 
beads and within bins of equal log width. The PSD of a bead 
in glycerol is shown for comparison. The indicated slopes 
are —1.75 and —2.0. Inset: Time series data for a 0.9 /im 
bead. For the smaller beads, deviations from scaling may 
indicate a failure of the continuum assumption. (B) : Real part 
(solid) and imaginary part (dotted) of the complex modulus 
G = G' (to) + iG"{ijj) calculated from the spectrum in (A) 
(5.0 /im bead). G' and G" scale as ^J^ with z = 0.76 ± 0.02 
determined from their ratio. The downturn in the spectra at 
about 20 kHz is due to an anti-aliasing filter. 

The fluctuation spectra of 0.9 /im silica beads in 2, 
2.5, and 3% crosslinked PAAm gels are shown in Fig. |^. 
The mesh size of this gel is about 50 A |2^ so that a 
continuum model easily applies. The 2% spectrum has 
a slope of about —1.5 near 100 Hz, consistent with the 
Rouse model and in clear contrast to the actin solutions. 
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The 2% spectrum decreases in slope below about 2 Hz 
due to the onset of an elastic plateau in G" (the other 
spectra reflect a plateau over most of their range, aside 
from the low-frequency instrumental noise from drift be- 
low a few Hz). The calculated G'{lu) are shown explicitly 
in Fig. 1^ lower inset, where the plateau values of G' arc 
2.0 Pa (2%), 24 Pa (2.5%) and 100 Pa (3%). Because the 
plateau extends to such high frequencies, G'{uj) does not 
yet exhibit the G{uj) ~ cj^/^ scaling Q as suggested by 
the PSD slope of —1.5. Close to the gelation threshold, 
the shear modulus increases very rapidly with concentra- 
tion, as also found by others . 
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FIG. 2. Power spectra of thermal motion of 0.9 fim beads 
in polyacrylamide gels of different concentrations, averaged as 
in Fig. cut off by an anti-aliasing filter at 20 kHz. Lower 
inset: The elastic modulus G'{uj) obtained from each of these 
spectra, as in Fig. |i]b. Upper inset: Plateau moduli (values 
at 1 Hz), as a function of concentration, of the three curves 
in the lower inset. 

The present technique differs in its spatial resolution 
and wide frequency range from prior microrheological 
studies such as those of Refs. |^-^ employing video mi- 
croscopy. The extended frequency range has allowed 
us to determine the high-frequency scaling of the shear 
modulus for actin well above the elastic plateau. We 
have described the inherent limitations in this and simi- 
lar microrheology studies at low frequencies. For rather 
sparse networks, the low-frequency response cannot be 
described by a shear modulus alone, but must also take 
into account the local compression of the network; this 
may invalidate the generalized Stokes formula, Eq. (|^), 
for frequencies accessible to video. In comparison to most 
commercial rheometers, our technique allows the study 
of soft materials and small sample volumes, which may 
be particularly useful in biological applications. 
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